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Denisonia devisi ACN-Den-1 [JAA74736] 
Echiopsis curia ACN-Ech-1 [JAA74766] 
Acanthophis wellsi ACN-Aca-1 [JAA74642] 
■ Suta fasciata ACN-Sut-1 [JAA74986I 
Ophiophagus hannah ache (VG) 
Bungarus fasciatus Ache [Q92035] 
Bungarus fasciatus Ache [AAC59905] 
Thamnophis elegans ache (Liver) 
[r Pantherophis guttatus ache transcript 1 (SAL, SCG, SK, Brain, VMNO) 
®i 3t Opheodrys aestivus ache transcript 2 (SCG) 
I utjDpheodrys aestivus ache transcript 1 (SAL, SCG) 
2£ U Echis coloratus ache transcript 1 (VG, SCG) 
b! Echis coloratus ache transcript 2 (SCG, SK) 
J Python regius ache transcript 2 (SCG) 
9^ Python regius ache transcript 1 (SAL, SCG, SK) 
' — Eublepharis macuiarius ache transcript 1 (SAL, SCG, SK) 
Gallusgallus Ache [AAA60456] 
Monode/pfi/sdomesf/caAche [XP 001371375] 
Homo sapiens Ache [AAA68151] 
Hartusnorveg/cus Ache [AAH94521] 
Mus musculus Ache [NP 033729] 
Xenopus tropicalis Ache [ENSXETG000000 17226] 
Damo rer/o Ache [ENSDARG00000031796] 
Takifugu rubripes Ache [ENSTRUG00000005195] 
Torpedo californica Ache [CAA27169] 
Lepisosteus oculatus Bche [ENSLOCG00000003535] 

Xenopus fropjca/j'sCholinesterase-like [XP 002931610] 
1 00r ftfas musculus Bche rENSMUSG000000277921 
1 02 L fiatfi/s norvegicus Bche [ENSRNOG00000009826] 
^ Homo sapiens Bche [ENSG000001 14200] 

' Monodelphis domestica Bche [ENSMODG00000016264] 

Ec/ijs co/orarus bche (VG, SCG, SK) 
Python regius bche (SCG) 
Anolis carolinensis bche [ENSACAG0000000771 9] 




ache 
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■ Eublepharis macuiarius bche (SAL, SCG) 
- Pelodiscus sinensis Bche [ENSPSIG00000007392] 

R1[ | Taeniopygia guttata Bche [ENSTGUG00000011088] 

1 0Oll — Columba livia Cholinesterase isoform X1 [XP 00551 301 2] 
83l_T Gallusgallus Bche [ENSGALG00000009478] 
life Meleagrisgallopavo Bche [ENSMGAG00000010452] 
1 0O r Anolis carolinensis cholinesterase-like [XP 003217998] 
l-Ano//scaro/;nens/s cholinesterase -like [XP 003217999] 

Eublepharis macuiarius cholinesterase-like (SK) 
Eublepharis macuiarius cholinesterase-like (SCG) 
Pantherophis guttatus cholinesterase-like (SCG, SK) 

Thamnophis elegans cholinesterase-like (Pooled tissue) 
Ophiophagus hannah Ache [ETE57261] 
Ovophis okinavensis Ache [BAN89435] 
Python molurus cholinesterase-like (Pooled tissue) 
100 i — Ophiophagus hannah Ache [ETE57655] 

' Thamnophis elegans c holi nesterase-li ke (Liver, Pooled tissue) 

- Anolis carolinensis Cholinesterase-like [XP 003218003] 
■ Python regius cholinesterase-like (SAL, SCG, SK) 
- Echis coloratus cholinesterase-like (VG, SCG) 
1 00LT Ophiophagus hannah Cholinesterase [ETE70564] 
gg|r Opheodrys aestivus cholinesterase-like (SAL, SCG, SK) 
1 00 1 - Pantherophis guttatus cholinesterase-like (SAL, SCG, SK) 
- Latimeria chalumnaeCEL [EN SLACG000000 15289] 
100 i — Monodelphis domestica CEL [ENSMODG00000012628] 
' — Homo sapiens CEL [ENSG000001 70835] 

Anolis carolinensis CEL [ENSACAG00000004871] 
Pelodiscus sinensis CEL [ENSPSIG00000012415] 
Taeniopygia guttata CEL [ENSTGUG00000005306] 
Gallus gallus CEL [ENSGALG00000003471] 
Meleagrisgallopavo CEL [ENSMGAG00000006148] 




che-like 
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Supplementary figure SI 

Maximum likelihood tree of acetylcholinesterase (ache), butyrylcholinesterase (bche) and cholinesterase-like sequences. 
Genbank accession numbers or Ensembl gene IDs are given in square brackets. Sequences from this study are in orange text 
and the location of transcripts is provided in bracketed blue text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, 
scent gland; VMNO, vomeronasal organ). Numbers above branches correspond to Bootstrap values for 500 replicates. 
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521- Micropechis ikaheka Cystatin E/M [E3P6P3] 
99^ Pseudechis australis [239618542 F] 
& A\Pseudonaja textilis [239618534 F J] 

ihOxyuranus scutellatus Cystatin [ACR83842] 
S -Ophiophagus hannah cystatin -e/m (VG, AG) 
Waya atra Cystatin [P81714] 
Naja kaouthia Cystatin [ACR83850] 

h Thamnophis efegans cystatin-e/m (Liver, Pooled tissue) 

94 j^r- Opheodrys aestivus cystatin-e/m (SAL, SCG, SK) 

7?— Pantherophis guttatus cystatin-e/m (SAL, SCG, SK, VMNO) 

r—Crotalus adamanteus Cystatin [AFJ49804] 
rdj— Echis coloratus cystatin-e/m (VG, SCG, SK) 
72LJ— Bitis gabonica Cystatin [AY430403] 
8? — Bitis arietans Cystatin [P08935] 

57 r Pantherophis guttatus cystatin (Brain) 

91I ' Ophiophagus hannah cystatin (AG) 

Pantherophis guttatus cystatin (SAL, SCG, SK) 

36? Opheodrys aestivus cystatin (SAL, SCG, SK) 

Thamnophis eiegans cystatin (Pooled tissue) 
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■ Thamnophis eiegans cystatin (Pooled tissue) 
Thamnophis eiegans cystatin (Pooled tissue) 



Thamnophis eiegans cystatin (Pooled tissue) 

Thamnophis eiegans cystatin (Pooled tissue) 
Pantherophis guttatus cystatin (SCG) 
Opheodrys aestivus cystatin (SCG) 

— Ophiophagus hannah cystatin (VG) 
Python regius cystatin e/m (SAL, SCG, SK) 
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Pogona barbata Cystatin [DQ1 39910] 
-Iguana iguana Cystatin [JX467151] 
Anolis carolinensis cystatin e/m [ENSACAG00000014558] 

Eublephans macularius cystatin e/m (SAL, SCG, SK) 
Chameleo catyptratus Cystatin E/M [JX467152] 



r-Homo sapiens Cystatin E/M [ENSG00000175315] 
LJOOr Mus musculus Cystatin E/M [ENSMUSG00000024846] 
L Rattus norvegicus Cystatin E/M [ENSRNOG00000020455] 
• Monodeiphis domestica Cystatin E/M [ENSMODG00000028390] 
Latimeria chalumnae Cystatin E/M [ENSLACG00000022489] 
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Supplementary figure S2 

Maximum likelihood tree of cystatin-e/m sequences. Genbank accession numbers or Ensembl gene IDs are given in square 
brackets. Sequences from this study are in orange text and the location of transcripts is provided in bracketed blue text (VG, 
venom gland; SAL, salivary gland; SK, skin; SCG, scent gland; AG, accessory gland; VMNO, vomeronasal organ). 
Numbers above branches correspond to Bootstrap values for 500 replicates. 
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Opheodrys aestivus cystatin-f (SAL, SCG, SK) 
Pantherophis guttatus cystatin-f (SAL, SCG, SK) 



99 



— Echis coloratus cystatin-f (VG, SCG) 
1 — Python regius cystatin-f (SAL, SCG, SK) 
■ Eublepharis macularius cystatin-f (SAL, SK) 
Gallus gallus Cystatin F [ENSGALG00000008660] 



Taeniopygia guttata Cystatin F [ENSTGUG00000006Q40] 

■ Monodelphis domestica Cystatin F [ENSMODG00000006300] 



100 



Homo sapiens Cystatin F [ENSG00000077984] 

I — Rattus norvegicus Cystatin F [ENSRNOG00000006767] 



1001 — Musmusculus Cystatin F [ENSMUSG00000068129] 
■Xenopus tropicalis Cystatin F [ENSXETG00000012775] 



■Latimeria chalumnae Cystatin F [ENSLACG00000002690] 
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Supplementary figure S3 

Maximum likelihood tree of cystatin-f sequences. Genbank accession numbers or Ensembl gene IDs are given in square 
brackets. Sequences from this study are in orange text and the location of transcripts is provided in bracketed blue text (VG, 
venom gland; SAL, salivary gland; SK, skin; SCG, scent gland). Numbers above branches correspond to Bootstrap values 
for 500 replicates. 
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— Ophiophagus hannah esp-e1 (VG, AG, PT) 
Thamnophis elegans esp-e1 (Liver) 
Opheodtys aestivus esp-e1 (SAL, SCG, SK) 
Pantherophis guttatus esp-e1 (SAL, SCG, SK, Brain, VMNO) 
Echis coloratus esp-e1 (VG, SCG, SK) 

LCrotalus adamanteus ESP E1 [AFJ49997] 
i 
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■Crotalus horridus epididymal secretory protein E1-like protein [JAA97132] 
-Python regius esp-e1 (SAL, SCG, SK) 
■ Python molurus esp-e1 (Pooled tissue) 

Eublepharis macularius esp-e1 (SAL, SCG, SK) 

Anolis carolinensis esp-e1 [XP 003214431] 

■Pogona vitticeps esp-e1 (Brain) 
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Pelodiscus sinensis epididymal secretory protein E1 [XP 006117809] 
Taeniopygia guttata ESPE1 [ENSTGUG00000016271] 
Gallus gallus ESP E1 [NP 001026374] 



TP- Meleagris gallopavo ESP E1 [ENSMGAG00000012679] 
Womo sapiens ESP E1 [NP 006423] 
Mus musculus ESP [BAA35183] 
Rattus norvegicus ESP E1 [NP 775141] 



■Latimeria chalumnae ESP E1 [ENSLACG00000018729] 



Supplementary figure S4a 

Maximum likelihood tree of epididymal secretory protein el (esp-el) sequences. Genbank accession numbers or Ensembl 
gene IDs are given in square brackets. Sequences from this study are in orange text and the location of transcripts is 
provided in bracketed blue text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, scent gland; AG, accessory gland; 
PT, pooled tissue; VMNO, vomeronasal organ). Numbers above branches correspond to Bootstrap values for 500 replicates. 
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Ornithorhynchus anatinus [ENSOANG00000030580] 
Chrysemys picta bellii epididymal sperm-binding protein 1 -like [XP 005311801] 
Ornithorhynchus anatinus [ENSOANG00000011916] 
Ornithorhynchus anatinus [ENSOANG00000011918] 
•Liophis poecilogyrus vMMP-Liol [ABU68569] 

Ophiophagus hannah Epididymal sperm-binding protein 1 [ETE57306] 
Varanus komodoensis ESP-Vkom1 [ADK39285] 
Varanus komodoensis matrix metalloprotease type 2 Var2 [ABY89671] 
Varanus komodoensis ESP-Vkom2 [ADK39286] 
— Varanus komodoensis matrix metalloprotease type 2 Var1 [ABY89663] 
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Varanus gouldii ESP-Vgou1 [ADK39283] 
Varanus indicus ESP-Vind2 [ADK39287] 
^Anolis carolinensis [ENSACAG00000029230] 
Eublephairs macularius (SCG) 
rCylindrophis ruffus ESP-Cyl-1 [AGI97178] 
|i/4/JO//s carolinensis [ENSACAG00000029124] 
To&Anolis carolinensis (Pooled tissue) 
^Anolis caroiinenesis [ENSACAG00000028839] 

^r/4no//s carolinensis epididymal sperm-binding protein 1 -like [XP 003224568] 
^Anolis carolinensis [ENSACAG00000010748] 

Ornithorhynchus anatinus [ENSOANG00000000388] 



Anolis carolinensis [ENSACAG00000010744] 
Supplementary figure S4b 

Maximum likelihood tree of epididymal secretory protein sequences. Genbank accession numbers or Ensembl gene IDs are 
given in square brackets. Sequences from this study are in orange text and the location of transcripts is provided in 
bracketed blue text (SCG, scent gland). Numbers above branches correspond to Bootstrap values for 500 replicates. 
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93 p. Opheodrys aesfrVus ficolin (SAL, SCG) 

Pantherophis guttatusf\co\\n-a (SAL, SCG, SK) 
- Pantherophis guttatusfico\\n-b (SAL, SK) 
— Philodryas olfersii Ficoli n[AFH66977] 
Thamnophis elegans ficolin (Pooled tissue) 
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'Cerberus rynchops Ryncolin-4 [D8VNT0] 
1 Cerberus rynchops Ryncolin-1 [D8VNS7] 
■ Cerberus rynchops Ryncolin-2 [D8VNS8] 
■ Cerberus rynchops Ryncolin-3 [D8VNS9] 
■ Enhydris polylepis Ficolin [AFH66974] 
•Echis coloratus ficolin (VG, SCG, SK) 



ioo I Leroe ' 

^ 1 99 f Cerfaer 

wCerb 

69l L_, 



i— Cylindrophis ruffus Ficolin [JX467168] 
30l9a— Morelia spilota Ficolin [JX467167] 

Python regius ficolin (SAL) 
gg ^ Opheodrys aestivus ficolin (SCG) 

i— TTiamnop/i/se/ega/isficolin (Pooled tissue) 
7£ Mi Opheodrys aestivus ficolin (SAL) 

9 *l Or 
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Opheodrys aestivus ficolin (SK) 
— Varanus fromodoens/sVeficolin [ADK46899] 
■ Eublepharis macularius ficolin (SAL, SCG, SK) 
Pogona vitticeps ficolin (Brain) 
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' Chlamydosaurus fr/ng/'/Veficolin [AGI97188] 
•Anolis carolinensis ficolin [ENSACAG00000025687] 
•Anolis carolinensis ficolin [ENSACAG00000017495] 
100 ^— Mus musculus Ficolin 1 [NP032021] 
■— Rattus norvegicus Ficolin 1 [NP1 12638] 
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Homo sapiens ficolin 1[NP001994] 
Chelonia mydas Ficolin-1 [EMP29986] 

— Columbalivia Ficolin-1 [EMC88538] 

99 Taeniopygia guttata ficolin-1 -like [XP002189999] 

— Galhisgallus Ficolin-1 (Brain) 

— ^Afias platyrhynchos Ficolin-1 [EOB04912] 
— Laf/meriacfrafarnnaeficolin-like [ENSLACG00000015719] 
■ Latimeria cha/umnae ficolin-like [ENSLACG00000015622] 



• Xenopus tropicalisfcnl [NP001 025596] 



• Xenopus tropicalis TNR [ENSXETG00000009195] 
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Supplementary figure S5 

Maximum likelihood tree of ficolin sequences. Genbank accession numbers or Ensembl gene IDs are given in square 
brackets. Sequences from this study are in orange text and the location of transcripts is provided in bracketed blue text (VG, 
venom gland; SAL, salivary gland; SK, skin; SCG, scent gland). Numbers above branches correspond to Bootstrap values 
for 500 replicates. 
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Opheodrys aestivus klk (SCG) 
Pantherophis guttatus klk (SCG) 
Echis coloratus klk-a (VG, SCG, SK) 
Python regius klk (SCG, SK) 

- Eublepharis macularius klk (SAL, SCG) 
Monodelphis domestica KLK 15 [ENSMODG00000006460] 

- Homo sapiens KLK1 5 [ENSG00000174562] 
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Rartt/s norvegicus KLK15 [ENSRNOG00000033805] 
JWt/s musculus KLK1 5 [ENSMUSG00000055193] 
Monodelphis domestica KLK1 1 [ENSMODG00000006536] 
■ Homo sapiens KLK1 1 [ENSG000001 67757] 
Mus musculus KLK1 1 [NP 001170844] 
Rartt/s norvegicus KLK1 1 [EN SRNOG000000 18742] 
100 |— Mus musculusKL.Ki4[NP 777355] 

100 



Rartt/s norvegicus KLK1 4 [ENSRNOG00000033706] 
Homo sapiens KLK14 [EN SG00000 129437] 



. Monodelphis domestica KLK14 [ENSMODG00000020815] 
- Ecm's coloratus klk-b (SCG, SK) 

Thamnophis elegans (Pooled tissue) 
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Homo sapiens KLK1 [AAH05313] 
Mus musculus KLK 1 [AAA3934] 



Rattus norvegicus KLK1 [NP1 1371 1] 



■ Human KLK8 [ENSG000001 29455] 



99 j— Mouse KLK8 [ENSMUSG00000064023] 



■ Rat KLK8 [ENSRNOG0000001 8580] 
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Supplementary figure S6 

Maximum likelihood tree of kallikrein (klk) sequences. Genbank accession numbers or Ensembl gene IDs are given in 
square brackets. Sequences from this study are in orange text and the location of transcripts is provided in bracketed blue 
text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, scent gland). Numbers above branches correspond to 
Bootstrap values for 500 replicates. 
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■Latimeria chalumnae Kunitzl [ENSLACG0000001 5634] 
r—Xenopus tropicalis Kunitzl [ENSXETG00000021555] 

iMus musculus Kunitzl [ENSMUSG00000027315] 
^Rattus norvegicus Kunitzl [ENSRNOG00000012811] 
jflHomo sapiens Kunitzl [ENSG00000166145] 
'■Monodelphisdomestica Kunitzl [ENSMODG00000000212] 
Meleagris gallopavo Kunitzl [ENSMGAG0000001 1 126] 
Gallusgallus Kunitzl [ENSGALG00000008469] 
Taeniopygia guttata unknown [ENSTGUG00000010764] 

■Anolis carolinensis kunitzl [ENSACAG00000003408] 
Python regius kunitzl (SAL, SCG, SK) 
Python regius kunitzl (SCG) 
Ophiophagus hannah kunitzl (AG, PT) 
Thamnophis elegans kunitzl (Liver) 
■Echis coloratus kunitzl (VG, SCG, SK) 
Opheodrys aestivus kunitzl (SAL, SCG, SK) 
Pantherophis guttatus kunitzl (SAL, SCG, SK) 
Opheodrys aestivus kunitz2 (SAL, SCG, SK) 
Pantherophis guttatus kunitz2 (SAL, SCG, SK) 

Austrelaps labialis Putative Kunitz-type proteinase inhibitor [B2BS84] 
Ophiophagus hannah kunitz2 (VG, AG, PT) 
Micrurus fulvius Kunitz inhibitor 4 [JAB52835] 
s^Thamnophis elegans kunitz2 (Liver) 
" is coloratus kunitz2 (VG, SCG, SK) 
ion regius kunitz2 (SAL, SCG, SK) 
50|rHefoderma suspectum cinctum kunitz-Hel-1 [AGI97151] 
^Abronia graminea kunitz-Abr-1 [AFU63209] 
•Iguana iguana kunitz-lgu-1 [AGI97181] 
-Chamaeleo caiyptratus kunitz -Cha-1 [AGI971 79] 
t Anolis carolinensis kunitz2 [H9GQ81] 
T|Afio/js carolinensis kunitz2 [H9G8B3] 
Va/jo/j's carolinensis kunitz2 [H9GRC6] 
Eublepharis macularius kunitz2 (SAL, SCG, SK) 
Chelonia mydas Putative Kunitz-type proteinase inhibitor [EMP28254] 
Monodelphis domestica Kunitz2 [ENSMODG00000013117] 
■Homo sapiens Kunitz2 [ENSG00000167642] 
JWusmuscu/usKunitz2[ENSMUSG00000074227] 
Rattus nofveg/ci/sKunitz2 [ENSRNOG00000020636] 
-Latimeria chalumnae Kunitz [ENSLACG00000003368] 
\Homo sapiens Ku nitz3 [ENSG000001 01446] 
Gorilla gorilla kunitz-type protease inhibitor 3 [XP 004062304] 
Callithrix jacchus kunitz-type proteaseinhibitor [3 XP 002747640] 
Musmi/sco/usKunitz3 [ENSMUSG00000074596] 

•Anolis carolinensis ku nitz 3 [ENSACAG00000023100] 




S2r-Rattus norvegicus Kunitz4 [ENSRNOG00000014794] 
SSP-Mus musculus Ku nitz4 [ENSMUSG0000001 731 0] 

Homo sapiens Kunitz4 [ENSG00000149651] 
X\Mus musculus Kunitz type 5 precursor [IMP 001035144] 

Rattus norvegicus Kunitz type 5 precursor [NP 001 008881] 
vxtPhilodryas olfersii Kunitoxin-PhM [ABU68488] 
^^Pantherophis guttatus kunitz (SAL, SCG, SK) 

ll^Chamaeleo calvotratus kunitz-Cha-2 [AGI97180] 
W'— Python regius kunitz (SCG) 
^Eublepharis macularius kunitz (SAL, SCG, SK) 
TAustrelaps superbus su perbin-4 [B5KL41 ] 
rVermicellaannulata KP-Ver-2 [JAA74997] 
P7ropf'dec/]/scarinaf«scarinatin-2 [B5KL33] 

L Sufa nigriceps Kunitz-type serine protease inhibitor isoform 1 [ACY68697] 
■(Va/fer/n/ies/aaegyptfaKunitzinhibitor I [ABX82867] 
f Micrurus fulvius Kunitz inhibitor 6 [JAB52833] 
^Micrurus fulvius Kunitz inhibitor 8 [JAB52831] 
r Micrurus fulvius Kunitz inhibitor 7 [JAB52832] 
^Walterinnesia aegyptia Kunitz inhibitor IV [ABX82870] 
■Vipera nikolskii Kunitz type inhibitor [CBW30779] 
-Micrurus fulvius Kunitz inhibitor 3 [JAB52836] 
rBitis gabonica Kunitz protease inhibitor 2 [AAR24535] 
^Bitis gabonica Kunitz protease inhibitor 1 [AAR24526] 
lfli/srretepssj;peri)asKunitz/BPTI-like toxin [B5L5M7] 
'Telescopusdhara Kunitoxin-Tel1 [ABU68487] 
'Oxyuranus sc«fe//afi/staicatoxin [ACC77765] 
Opheodrys aestivus kunitz (SCG, SK) 
Pantherophis guttatus kunitz (SK) 
-Opheodrys aestivus kunitz (SCG) 
-Echis coloratus kunitz (VG, SCG, SK) 
~Dendroaspis polylepis Dendrotoxin-1 [P00979] 
^Dendroasp/spo/ylep/sdendrotoxin K [AAB26998] 
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Supplementary figure S7 

Maximum likelihood tree of kunitz 1 and 2 sequences. Genbank accession numbers or Ensembl gene IDs are given in square 
brackets. The three kunitz2 sequences from the Green Anole (Anolis carolinensis) are derived from the genome of this 
species as used previously by Fry et al. [2]. Sequences from this study are in orange text and the location of transcripts is 
provided in bracketed blue text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, scent gland; AG, accessory gland; 
PT, pooled tissue). Numbers above branches correspond to Bootstrap values for 500 replicates. 
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Pantherophisguttatus\\pa-a (SAL) 
PantherophisguttatusUpa-a (SK) 
Pantherophisguttatus tipa-a (SCG) 
Opheodrysaestivus lipa-a (SAL, SCG, SK) 

Ec/i/sco/oratus lipa-a (VG, SCG, SK) 
Thamnophis elegans lipa-a (Pooled tissue) 
Python regius lipa-a (SAL, SCG, SK) 
Python molurus lipa-a (Pooled tissue) 
•Anoliscarolinensis lysosomal acid lipase [XP 003229744] 
^^^—Pogona vitticepsUpa (Brain) 
•Eublepharis macularius lipa-a (SAL, SCG, SK) 

Crocodylus niloticus lipa (Brain) 
Chelonia mydas Lysosomal acid lipase/cholesteryl ester hydrolase [EM P2461 5] 
Pelodiscussinensis LIPA [ENSPSIG00000016953] 
Coiumba livia Lysosomal acid lipase/cholesteryl ester hydrolase[EMC89691] 
Taeniopygia guttata lysosomal acid lipase/cholesteryl ester hydrolase [XP 004174885] 
Anas platyrhynchos | Lysosomal acid lipase/cholesteryl ester hydrolase [EOA94419] 
99 i~ Callus gallus lysosomal acid lipase/cholesteryl ester hydrolase [XP421661] 

Meieagrisgallopavo lysosomal acid lipase/cholesteryl ester hydrolase [XP003208075] 
83^^^^^^^— Gerrhonotus infernalis lysosomal acid lipase [ADK5561 0]| 
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•Anolis carotinensis\ysosoma\ acid lipase/cholesteryl ester hydrolase [XP003224429] 

•PantherophisguttatusWpa (SCG) 
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■ Varanusg/auerr/lysosomal acid lipase [ADK55609] 
Varanusscataris lysosomal acid lipase [ADK55608] 

■ Varanuskomodoensis lysosomal acid lipase [ADK55607] 

97 ■ Eublepharis maculariustipa (SAL) 

Eublepharis macularius lipa (SCG, SK) 
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— Anolis carotf/Tens/slysosomal acid lipase/cholesteryl ester hydrolase [XP003224425] 
■ Pantherophisguttatus lipa (SCG) 
— Eublepharis macularius lipa (SCG, SK) 
96 1 ^— ^— OpfteodVysaeslrais lipa (SCG) 
^—Pantherophisguttatus lip (SAL) 
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™Ec/7(sco/oraruslipa-b(VG, SCG) 
' Python regius\\pa-b (SCG) 
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•Anolis carolinensis\ysosoma\ acid lipase/cholesteryl ester hydrolase [XP 003224417] 
"Pylhonreg/uslipa (SCG) 
■ Eublepharis macularius lipa (SCG) 



^Eublepharis macularius lipa (SAL) 
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Eublepharis macularius lipa (SK) 
PantherophisguttatusWpa (SCG) 
~ Python regius lipa (SK) 

Monodelphis domestica LIPA [ENSMODG00000019646] 



93 ^— Homosapiens lysosomal acid lipase [AAB60328] 

' ,,^^99^^^^— Musmusculus lysosomal acid lipase/cholesteryl ester hydrolase [HP 067435] 
1 Rattus norvegicus lysosomal acid lipase/cholesteryl ester hydrolase [NP 036864] 
^^—Xenopus tropicalis lysosomal acid lipase/cholesteryl ester hydrolase [NP001015847] 
■ Latimeria chalumnae LIPA [ENSLACG00000000424] 
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Supplementary figure S8 

Maximum likelihood tree of lysosomal acid lipase (lipa) sequences. Genbank accession numbers or Ensembl gene IDs are 
given in square brackets. Sequences from this study are in orange text and the location of transcripts is provided in 
bracketed blue text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, scent gland). Numbers above branches 
correspond to Bootstrap values for 500 replicates. 
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—Bungarus mutticinctus NGF [AAB25729] 
rDrysdalia coronoides NGF [AEH95535] 
>— Denisonia devisi NGF-Den-2 [JAA74726] 
r Ophiophagus hannah ngf (VG) 
'Ophiophagus hannah ngf (Genomic- scaffoldl 2849.1) 
A/a/a sputatrix NGF [Q5YF89] 
Naja naja Venom nerve growth factor [P01140] 
Naja kaouthia nerve growth factor [A59218] 
Naja atra Venom nerve growth factor [P61898] 
^Walterinnesia aegyptia venom nerve growth factor [AEH59581] 
Micrurus fulvius NGF-1a (VG) 
Micrurus fulvius NGF-1b (VG) 
rNaja sputatrix NGF [AAS94268] 

Wa/a kaouthia nerve growth factor beta polypeptide [ACC85831] 
Ophiophagus hannah hypothetical protein L345 10532 [ETE63701] 
Ophiophagus hannah ngf (Genomic- scaffold2616.1) 
•Laticauda colubrina nerve growth factor beta polypeptide [ACC85824] 
Pseudonaja textilis NGF [Q3HXY9] 
-Oxyuranus microlepidotus NGF [ABA60117] 

ITrimeresurus flavoviridis NGF [BAG16511] 
Lachesis stenophrys nerve growth factor beta polypeptide [ACC85821] 
Ovophis okinavensis NGF [BAN82142] 
•Crotalus durissus NGF [Q9DEZ9] 
rBothrops jararacussu NGF [AAG12169] 
j^Bothrops asper nerve growth factor beta polypeptide [ACC85807] 
rMacrovipera lebetina NGF [P25428] 

J-Daboia russelii nerve growth factor beta polypeptide [ACC85813] 

Lec/i/s coloratus ngf (VG, SCG, SK) 
-Azemiops feae NGF-AZE1 [AAZ75635] 

^Thamnophis marcianus nerve growth factor beta polypeptide [ACC85837] 
jiOpheodrys aestivus ngf (SCG, SK) 
\Pantherophis guttatus ngf (SAL, SCG) 

'-Lampropeltis getula nerve growth factor beta polypeptide [ACC85822] 

Python regius ngf (SK) 

Python molurus ngf [AEQU01 0205372] (Genomic- Contig26560677) 
-Boa constrictor nerve growth factor beta polypeptide [ACC85806] 
r— Ano//s carolinensis ngf [ENSACAP00000003352] 
|r Varanus acanthurus NGF-VAR1 [AAZ75634] 

*Varanus salvator nerve growth factor beta polypeptide [ACC85794] 
^—Eublepharis macularius ngf (SAL, SCG, SK) 

ggi Taeniopygia guttata NGF [ENSTGUG00000000958] 

1 ipprGa/Zt/S gallus NGF [ENSGALG00000002556] 

meieagris gallopavo [NGFENSMGAG00000015674] 
Monodelphis domestica NGF [ENSMODG00000004863] 
m j-Homo sapiens NGF [ENSG000001 34259] 
lOorMus musculus NGF [AAA37686] 

Rattus norvegicus NGF [ENSRNOG00000016571] 
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•Latimeria chalumnae NGF [ENSLACG00000018594] 
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Supplementary figure S9 

Maximum likelihood tree of nerve growth factor (ngf) sequences. Genbank accession numbers or Ensembl gene IDs are 
given in square brackets. Sequences from this study are in orange text and the location of transcripts is provided in 
bracketed blue text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, scent gland). Numbers above branches 
correspond to Bootstrap values for 500 replicates. 
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■ Leioheterodon madagascariensis PLA 2 Group HE [AFH66961] 
Dispholidus typus PLA 2 Group HE [AFH66959] 



Pantherophis guttatus PLA 2 Group HE (SAL) 
Opheodrys aestivus PLA 2 Group HE (SAL) 

Echis coloratus PLA 2 Group HE (VG) 



98 p Python regius PLA 2 Group HE (SAL) 

I— Python regius PLA 2 Group HE (SAL, SCG) 
— Eublepharis macularius PLA 2 Group HE (SAL) 
•Pelodiscus sinensis PLA 2 Group HE [ENSPSIG00000005789] 

Taeniopygia guttataPLA 2 Group HE [ENSTGUG00000015886] 
Gallus gallus PLA 2 Group HE ENSGALG00000014176 

Meleagris gallopavo PLA 2 Group HE [ENSMGAG00000017176] 
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• Monodelphis domestica PLA 2 group HE [ENSMODG00000015910] 
'Mus musculus group HE secretory phospholipase A 2 [NP036174] 
• Homo sapiens PLA 2 group HE [ENSG000001 88784] 

^Xenopus laevis phospholipase A 2 group HE [NP001 087269] 



■ Danio rerio secretory phospholipase A 2 [NP001002350] 
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Supplementary figure S10 

Maximum likelihood tree of Group HE Phospholipase A 2 {PLA 2 Group HE) sequences. Genbank accession numbers or 
Ensembl gene IDs are given in square brackets. Sequences from this study are in orange text and the location of transcripts 
is provided in bracketed blue text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, scent gland). Numbers above 
branches correspond to Bootstrap values for 500 replicates. 
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m Monodetphis domestics ly-6/neurotoxin-like protein 1 -like [XM 001 381754] 
\—Mus muscutus Ly-GV neurotoxin homo log [AF141377] 
Python molurus 3ftx-a {Pooled tissue} 
Thamnophis elegans 3ftx-a (Pooled tissue) 
^ Pantherophis guttatus 3ftx-a (SAL) 
^Pantherophis guttatus 3ftx-a (SK) 
rGpheodtys aestivus 3ftx a (SK) 
Opheodrys aestivus 3ftx-a (SAL) 
^Pantherophis guttatus 3 fix a ( SCG) 
EcMs coJorafris 3ftx-a(VG, SCG, SK) 
'fytfjon regius 3ftx-a ( SAL) 

^-Lachesis muta three finger toxin-like precursor [ABD52883] 
Thamnophis elegans 3ftx (Pooled tissue) 
Thamnophis elegans 3ftx [Pooled tissue) 
; r7hrasops jacksoni 3FTx-Thr2 [EU 0296 84] 
H- Thrasops jacksoni 3FTx-Thr3 [E U 0 29 6 85] 
sa tThrasops jacksoni 3FTx-Thr1 [E U 0296 82] 
^Dispholidus typus 3FTx-Dis2 [EU029681] 
^Dispholidus typus 3FTx-Dis3 [EU 0296 8 3] 
" mr Boiga dendrophitia Denmotoxin [Q05ZW0] 
xp-Tetescopus dhara 3FTx-Tel4 [EU029686] 
log Tetescopus dhara 3FTx -Te 12 [E U 0 29 6 7 9] 
\Teiescopus dhara 3FTx-Tel3 [EU029G80] 
■Disphotidus typus 3FTx-Dis1 [EU029674] 
^fcej— Liophis poecilogyrus 3FTx Lio3 [ABUG8473] 
— i— Liophis poecilogyrus 3FTx-Lio2 [ABUG8472] 
r-Tetescopus dhara 3FTx Tell [EU029671] 
ILopheodrys aestivus 3ftx(SAL, SCG) 
iPython regius 3ftx (SCG) 
^Opheodrys aestivus 3ftx (SAL) 
T Opheodrys aestivus 3ftx (SCG) 
^Python regius 3ftx ( SCG) 
V Thrasops jacksoni 3FTx-Thr5 [EU036635] 
3 ^Disphotidus typus 3FTx-Dis4JEU03G63G] 
\-Trimorphodon biscutatus 3FTx-Tri1 [ABU G 847 5] 
-Leioheterodon madagascariensis 3FTx -Lei1 [E U 029 676] 
—Psammophis mossambicus 3FT x-P sal [EU0296G9] 

Sistrurus catenatus edwardsi three finger toxin 3 [DQ4G4283] 
ioa \Echis coioratus 3ftx-b (VG, SCG) 
S^Opheodrys aestivus 3ftx-b ( SCG) 

~f i- Sistrurus catenatus edwardsi three finger toxin 2 [DQ464282] 
^1 r Sistrurus catenatus edwardsi three finger toxin 1 [DG464281] 
^j-Sfsfrurt/s cafenafi/s edwardsi three finger toxin 4B [DQ464285] 
^S/stfwrus catenatus edwardsi three finger toxin 4A [DQ464284] 
—Pantherophis guttatus 3ftx (SCG) 
s rPantherophis guttatus 3ft x (SAL) 
l fiantf?erophfs guttatus 3ftx (SCG) 
■ Thamnophis elegans 3ftx (Pooled tissue) 

TYjamnopfiis elegans 3ftx (Pooled tissue) 

Thamnophis elegans 3ftx (Pooled tissue) 
- Thamnophis elegans 3ftx (Pooled tissue) 
— Opheodrys aestivus 3ftx(SCG) 

Pantherophis guttatus 3ftx (SCG) 

Enhyrfris polylepis 3FTx-Enh1 [EU 0296 68] 

-Python regius 3ftx ( SAL) 

»r Wa/a sputatrix wNTX-9 weak neurotoxin precursor [AF098924] 
pNaja atra Full Probable weak neurotoxin NNAM2 [P60814] 
Y-Najaatra Full Probable weak neurotoxin HHAM1 [Q9YGI2] 
^■Bungarus ftaviceps non-conventional three finger toxin isoform 7 [ADF50022] 
— Bungarus multicinctus cardiotoxin-like protein [AJ007741] 
*[r—Bungarusflaviceps short-chain three finger toxin isoform 3 [ADF50011] 
!p— Ophiophagus hannah Three-Finger Neurotoxin [3HH7 A] 
'—Lapemis hardwickiiLonq neurotoxin 2 [A3FM53] 
r-Bungarus multicinctus neurotoxin -like protein [AJ 006135] 
^%rBungarusflaviceps non-conventional three finger toxin isoform 3 [ADF50018] 
'■Bungarus candidus candoxin [AY142323] 
sop Taeniopygia guttata prostate stem cell antigen [XP0021 87004] 
W- Gall us gall us prostate stem cell antigen [XP418414] 
«l '—Columba tivia Prostate stem cell antigen partial [EMC78508] 
9 I mrMus muscutus Prostate stem cell antigen [HM 02821 G] 

^■Rattus norvegicus prostate stem cell antigen precursor [NP001165577] 

1 Ga//us£a//usLY6/PLAUR domain containing 2 [XM 001231772] 

-Meleagris galtopavo lymphocyte antigen 6E-like [XP 003205371] 
rjiGatlus gatlus lymphocyte antigen G complex [XM 418413] 
"irMeleagris gattopavo lymphocyte antigen GE-like[XM 003205322] 

Taeniopygia guttata ste m ce II antig en 2 I i ke [N M 001 2456 35] 
Chelonia mydas Lymphocyte antigen 6E [EMP26638] 
^eAnolis carolinensis lymphocyte antigen GE like [XM 003223921] 

— Homo sapiens lymphocyte antigen G complex locus E variant [BAD96794] 
Gattusgaltus lymphocyte antigen 6E precursor [NP990106] 
Taeniopygia guttata lymphocyte antigen 6E like [XP0021 87037] 

—Pantherophis guttatus lymphocyte antigen 6E-like (SAL, SK) 
—Anotis carotinensis lymphocyte antigen EE-like [XP003223972] 
— Xenopus tropicals prostate stem cell arrtigen-like[XPO02941544] 



Supplementary figure Sll 

Maximum likelihood tree of three finger toxin (Sftx) genes. Genbank accession numbers or Ensembl gene IDs are given in 
square brackets. Sequences from this study are in orange text and the location of transcripts is provided in bracketed blue 
text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, scent gland). Numbers above branches correspond to 
Bootstrap values for 500 replicates. 
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Hoplocephalus stephensiivesaryn isoform 1 [ABW74877] 
Demansia vesf/g/afa vespryn isoform 1 [ABW74870] 
Pseudonajatextifis ves pry n isoform 2 [ABW74865] 
Tropidechis carinatus ves pry n isoform 1 [ABW74874] 
Drysdalia coronoides Vesp23 [AEH95534] 
Cryptophis nigrescens vespryn isoform 1 [ABW74872] 
Cryptophis nigrescens vespryn isoform 2 [ABW74873] 
Notechis scutatus vespryn isoform 1 [ABW74866] 
Pseudechis austrafis vespryn [ABW74876] 
m Oxyuranus scutellatus vespryn isoform 1 [ABW74868] 

■ Drysdalia coronoides Vesp21 [AEH95532] 
Pseudonajatextifis vespryn isoform 1 [ABW74864] 

■ Drysdalia coronoides Vesp22 [AEH95533] 
Oxyuranus scutellatusvesoryn isoform 2 [ABW74869] 

mOphiophagus hannah Ohanin precursor [AAR07992] 

- Wa/a/raot/fri/aThaicobrin [P82885.] 
Pseudechis porphyriacusvespryn isoform 1 [ABW74875] 
Opheodrys aestivus vespryn (SCG, SK) 

Echis coloratus vespryn [KF1 14032] (genomic - contig 243235) 
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61 
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i^^Lachesis muta Ohanin-like protein [Q27J48] 
Crotalus adamanteusvesotyn [A E J 32004] 
Pantherophis guttatus vespryn (SAL, SCG) 
Pantherophis guttatus vespryn (genomic - contig 470301 ) 
Boa constrictor vespryn-Boa-1 [GI97190] 
99 i Python regius vespryn (SCG, SK) 

Python molurus vespryn [AEQU01 01 88142] (genomic - contig 26526219) 
Eublepharis macularius vespryn (SCG, SK) 
^^^^^^^^^ Anolis carolinensis ohanin-like [XP003218014] 
Pogona barbata vespryn-POGLH [Q2XXL4] 

Anolis carolinensis tripartite motif -containing protein 7-like [XP 003218019] 
^^^^^^^^^^^^^™Ano//scaro//>jens/s zinc finger protein RFP-like [XP 003218059] 




~Chelonia mycras Guanine nucleotide-binding protein subunit beta-2-like 1 [EMP25897] 
1 Chelonia mydas Tripartite motif -containing protein 39[EMP25888] 
97 . Meleagris gallopavo tripartite motif protein 39 [ADU03785] 



69, 



i Gall us gall us tripartite motif-containing 39 [N PO0 1006 196] 

^^^—Homo sapiens butyrophilin subfamily 3 member A1 [NP 001138480] 

Monode/ph/sdomesf/catripartite motif-containing protein 39-like [XP 001374079] 
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Supplementary figure S12 

Maximum likelihood tree of vespryn sequences. Genbank accession numbers are given in square brackets. Sequences from 
this study are in orange text and the location of transcripts is provided in bracketed blue text (VG, venom gland; SAL, 
salivary gland; SK, skin; SCG, scent gland). Numbers above branches correspond to Bootstrap values for 500 replicates. 
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ITropidechis carinatus Carwaprin-b [B5L5P0] 
■ Demansia vestigiata Veswaprin-c [B5L5P6] 
Notechis scutatus Notewaprin-a [B5G6H4] 
Notechis scutatus Notewaprin-b [B5G6H5] 
Demansia vestigiata Veswaprin-b [B5L5P5] 
• Cryptophis nigrescens Nigwaprin-a [B5G6H1] 
J* Pseudechis porphyriacus Porwaprin-b [B5L5N2] 
Cryptophis nigrescens Nigwaprin-b [B5L5N6] 

rHopiocephalus stephensii Stewaprin-a [B5G6H3] 
- Austrelaps superbus Supwaprin-a [B5KGY9] 
—Pantherophis guttatus waprin-like (SAL, SCG) 
Python regius waprin-like (SK) 
•Eubiepharis mactilarius waprin-like (SAL, SCG, SK) 
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100 
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— Thamnophis elegans waprin-like (Liver, Pooled tissue) 

— Ophiophagus Hannah waprin-like (VG) 

— Opheodrys aestivus waprin-like (SCG, SK) 

— Pantherophis guttatus waprin-like (SCG) 

— Philodryas otfersii Waprin-Phi3 [ABU68546] 

74 1 Ophiophagus hannah waprin-like (AG) 

Tj4i Pantherophis guttatus waprin-like (SAL, VMNO) 

Opheodrys aestivus waprin-like (SAL) 
' Python regius waprin (SCG) 



' Python regius waprin (SK) 



30 
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Ophiophagus hannah waprin (AG, PT) 

- Philodryas olfersii Waprin-Phi1 [A7X4K1] 

- Thamnophis eiegans waprin (Liver) 

- Opheodrys aestivus waprin (SAL, SCG, SK) 
■ Pantherophis guttatus waprin (SAL, VMNO) 

- Thamnophis eiegans waprin (Liver, Pooled tissue) 
Naja nigricoilis Nawaprin [P605E9] 

Ophiophagus hannah waprin (VG) 
Rhabdophis tigrinus Waprin-Rha1 [A7X4J4] 
— Pantherophis guttatus waprin (VMNO) 

Thrasops jacksoni Waprin-Thr1 [A7X4I7] 

Liophis poecilogyrus Waprin-Lio1 [A7X4L4] 

Philodryas olfersii Waprin-Phi2 [A7X4K7] 

r Crotalus adamanteus Waprin [AFJ49269] 
I Echis coloratus waprin (VG, SK) 

Ovophis okinavensis Waprin [BAN89446] 

Eublepharis macularius waprin (SAL, SCG, SK) 
Chrysemys picta bellii Waprin-Phi1-like [XP 005304910] 

100 r Alligator sinensis WAP four-disulfide core domain protein 2 [XP 006035848] 

Alligator mississippiensis WAP four-disulfide core domain protein 2-like [XP 006260197] 
Columba livia Waprin-Phi1-like [XP 005499578] 
Columba livia WAP four-disulfide core domain protein 2 [EMC83403] 
Gallus gallus WAP four-disulfide core domain protein 2-like [XP 003642517] 
Meleagris gallopavo WAP four-disulfide core domain protein 2-like [XP 003212514] 
97^— Rattus norveglcus Wfdc2 [ENSRNOG00000014739] 
Musmusculus Wfdc2 [ENSMUSG00000017723] 

Monodelphis domestica WAP four-disulfide core domain 2 [ENSMODG00000016177] 
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' Homo sapiens HE4 protein [CAA44369] 
' Pan troglodytes Wfdc2 [ENSPTRG00000013549] 
■ Xenopus tropicalis WAP four-disulfide core domain protein 2-like [XP 002935659] 
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Supplementary figure S13 

Maximum likelihood tree of waprin sequences. Genbank accession numbers or Ensembl gene IDs are given in square 
brackets. Sequences from this study are in orange text and the location of transcripts is provided in bracketed blue text (VG, 
venom gland; SAL, salivary gland; SK, skin; SCG, scent gland; AG, accessory gland; PT, pooled tissue; VMNO, 
vomeronasal organ). Numbers above branches correspond to Bootstrap values for 500 replicates. The two Echis coloratus 
sequences are likely to be alleles or the result of variation in the two individuals used to generate the transcriptomes. 
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Echis coloratus vegf-c (VG, SCG, SK) 
Crotalus horridus VEGF-C [JAA95033] 
Opheodrys aestivus vegf-c (SAL, SCG, SK) 
Pantherophis guttatus vegf-c (SAL, SCG, SK) 
Thamnophis efegans vegf-c (Liver) 
Micrurus fulvius VEGF-C [JAB52939] 
; ^Ophiophagus hannah vegf-c (VG) 

Python regius vegf-c (SAL, SCG, SK) 
^Python molurus vegf-c (Pooled tissue) 
-Anolis carolinensis vegf-c (Pooled tissue) 
§ -Eublepharis macularius vegf-c (SAL, SCG, SK) 
— Pogona vitticeps vegf-c (Brain) 
•Taeniopygia guttata VEGF-C [ENSTGUG00000006451] 
p tMeleagrls gallopavo VEGC-C [ENSMGAG00000007088] 
gBGeffus gallus VEGF-C [ENSGALG00000010847] 

Monodelphis domestica VEGF-C [ENSMODG00000003799] 
Homo sapiens VEGF-C [NP 005420] 
Rattus norvegicus VEGF-C [NP 446105] 
KKf Mus musculus VEGF-C [AAC52984] 
I— Latimeria chalumnae VEGF-C [ENSLACG00000002200] 
— Latimeria chalumnae VEGF-D [ENSLACG00000015806] 
Mus musculus VEGF-D [NP 034346] 
Rattus norvegicus VEGF-D [AAK96008] 
, - Homo sapiens VEGF-D [NP 004460] 
1 — J Gallus gallus VEGF-D [AAM12733] 

J 1 Taeniopygia guttata VEGF-D [XP 002197500] 
92|r Eublepharis macularius vegf-d (SAL, SCG, SK) 
gill - Python regius vegf-d (SAL, SCG, SK) 
68"— Anolis carolinensis vegf-d [ENSACAG0000001241 9] 
Macrovipera lebetina svVEGF [P82475] 
Daboia russellii VEGF-F [ACN22046] 
Vipers aspis VEGF [PS3942] 
Vipera ammodytes VEGF-F [ACN22045] 
~chis coloratus vegf-f (VG) 
Bills arietans Barietin [ACN2203S] 
Protobothrops mucrosquamatus TM-VEGF [AAS07632] 

Protobothrops flavoviridis vascular endothelial growth factor-like protein [BANS2012] 
Ovophis okinavensis vascular endothelial growth factor-like protein [BAN82145] 
Agkistrodon piscivorus apiscin precursor [ACN22039] 
Bolhrops jararaca VEGF-F [AAK52103] 
Bothrops insularls VEGF-F [Q90X24] 
Crotalus adamanteus vegf-f (VG) 
Crotalus atrox cratrin precursor [ACN22040] 
Sistrurus catenatus edwardsi VEGF-F [BOVXV3] 
Mus musculus VEGF-B [AAB06273] 

Homo sapiens VEGF-B [NP 003368] 
Monodelphis domestica VEGF B-like [XP 001378447] 
Eublepharis macularius vegf-b (SAL, SCG, SK) 
Echis coloratus vegf-b (VG, SCG, SK) 
321 Python regius vegf-b (SAL, SCG, SK) 
* Python molurus vegf-b (Pooled tissue) 
Ophiophagus hannah vegf-b (AG, PT) 
Thamnophis eiegans vegf-b (Liver) 
Opheodrys aestivus vegf-b (SAL, SCG, SK) 
Pantherophis guttatus vegf-b (SAL, SCG, SK) 
95r-Homo sapiens VEGF-A [NP 001020537] 

Mus musculus VEGF-A [ENSMUSG00000023951] 
9a- Taeniopygia guttata VEGF-A [ENSTGUG00000008137] 
~ Gallus gallus VEGF-A [ENSGALG00000010290] 

45|t — Anolis carolinensis vegf-a [ENSACAG0000000161 2] 
Thamnophis eiegans vegf-a (Liver) 
Opheodrys aestivus vegf-a (SAL, SCG, SK) 
Python regius vegf-a (SAL, SCG SK) 
Python molurus vegf-a (Pooled tissue) 
Pantherophis guttatus vegf-a (SAL, SCG, VMNO) 
jl Pseudonaja modesta VEGF-Pse-3 [JAA74898] 

V Suta nigriceps VEGF-A [ACY68720] 

f Ophiophagus hannah vegf-a (VG) 

t Micrurus futvlus vascular endothelial growth factor 2 [JAB52756] 
I Micrurus fulvius vascular endothelial growth factor 3 [JAB52755] 
■ Micrurus fulvius vascular endothelial growth factor 1 [JAB52757] 

lEchis coloratus vegf-a (VG, SCG, SK) 
-Bitis gabonica VEGF-A [AAR06855] 
Crotalus adamanteus VEGF 1a [AFJ49265] 

Crotalus adamanteus vascular endothelial growth factor 1b [AFJ49266] 
t protobothrops flavoviridis vascular endothelial growth factor-like protein [BAN89382] 
Protobothrops flavoviridis vascular endothelial growth factor-like protein [BAN89381] 
Ovophis okinavensis vascular endothelial growth factor-like protein [BAN89443] 
Ovophis okinavensis vascular endothelial growth factor-like protein [BAN89442] 
Ovophis okinavensis vascular endothelial growth factor-like protein [BAN89445] 
Ovophis okinavensis vascular endothelial growth factor-like protein [BAN89444] 
•Latimeria chalumnae PGF [ENSLACG00000018931] 
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Supplementary figure S14 

Maximum likelihood tree of vascular endothelial growth factor (vegf) sequences. Monophyletic clades representing the five 
members of the vegf family (a-f) are indicated and the putative viper venom-specific ve gf-f clade is shaded green. Genbank 
accession numbers or Ensembl gene IDs are given in square brackets. Sequences from this study are in orange text and the 
location of transcripts is provided in bracketed blue text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, scent 
gland; AG, accessory gland; PT, pooled tissue; VMNO, vomeronasal organ). Numbers above branches correspond to 
Bootstrap values for 500 replicates. 
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100j Crotalus adamanteus LAAO [AFJ49193] 

9 n Crotalus atrox LAA0 [AAD45200] 
1 0Q V- Lachesis muta LAAO [AFP89360] 

4— Calloselasma rhodostoma LAAO [CAB71136] 
100 *— Viridovipera stejnegeri LAAO [AAQ16182] 
Daboia russellii siamensis LAAO [AAZ08620] 
- Echis ocellatus LAAO [CAQ72894] 
100i Echis coloratus laao-b1 (VG) 
99* Echis coloratus laao-b2 (VG) 
1 DOT" Opheodrys aestivus laao-b (SCG) 
Pantherophis guttatus laao-b (SCG) 
Naja atra LAAO [ABN72546] 
Oxyuranus scutellatus LAAO [AAY89680] 
Bungarus fasciatus LAAO [ABN72540] 
IQO r Python regius laao-b (SCG) 

* Python molurus [AEQUM 0327398] (Genomic- contig 26804723) 
^^—Anolis carolinensis l-ami no-acid oxidase-like isoform 1 [XP003217921] 




100 




■ Xenopus tropicalis l-amino-acid oxidase-like [XP002938651] 
— Xenopus tropicalis LAAO [XP002937668] 
Meleagris gallopavo LAAO [ADU03782] 
Gallus gallus LAAO [NP001092821] 

Taeniopygia guttata LAAO [XP0021 89436] 
Columba livia LAAO [EMC90054] 
Chelonia mydas LAAO [EMP41398] 
Ecftis coloratus laao-a (VG, SCG, SK) 
Ophiophagus hannah LAAO [ABN72538] 
Anotis carolinensis l-amino-acid oxidase-like [XP003217939] 



Eublepharis macularius laao (SAL, SCG) 

Eublepharis macularius laao (SAL, SCG, SK) 
Anolis carolinensis [ENSACAG00000028748] 
Physignathus lesueurii LAAO 
Homo sapiens D-amino-acid oxidase [NP001908] 



0.5 



Supplementary figure S15 

Maximum likelihood tree of l-amino acid oxidase (laao) sequences. Genbank accession numbers or Ensembl gene IDs are 
given in square brackets (the Australian Water Dragon {Physignathus lesueurii) sequence does not have a Genbank 
accession number and is derived from a sequence in [2]). Sequences from this study are in orange text and the location of 
transcripts is provided in bracketed blue text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, scent gland). 
Numbers above branches correspond to Bootstrap values for 500 replicates. 
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*r Ophiophagus hannah Opharin [ACN93671] (Thailand) 
Y Ophiophagus hannah crisp-a (Indonesia) 



i 



Bungarus candidus Bc-CRPb [ACE73577] 
Bungarus candidus Bc-CRPa [ACE73578] 
Naja atra natrin-2 [Q7ZZN8] 
Naja kaouthia kaouthin-2 [ACH73168] 
|- Oxyuranus microlepidotus CRISP-OXY1 [AAZ75602] 
Hydrophis hardwickii CRISP 2 [Q8UW11] 
Hydrophis hardwickii CRISP 1 [Q8UW25] 
ei|_r Pseudechis porphyriacus pseudecin [Q8AVA3] 
w 1 - Pseudechis australis pseudechetoxin [Q8AVA4] 



'-Laticauda semifasciata Latisemin [Q8JI38] 

j- Opheodrys aestivus crisp-a ( 
- fell Pantherophis guttatus crisp-a (\ 
9 ™ Pantherophis guttatus crisp-a (£ 
r-Echis coloratus crisp-a (VC 
y_i — Trimeresurus gracilis Og-CRPb [ACE73570] 
73 — Trimeresurus flavoviridis Serotriflin [P0CB15] 



- Boa constrictor crisp-a 

Agkistrodon piscivorus piscivorin [AA062994] 
Crotalus atrox catrin [AA062995] 
Sistrurus catenatus CRISP [ABG26992] 



(E00001412) 



h — Gloydius blomhoffi ablomin [AAM45664] 
s<J— Protobothrops jerdonii CRISP [AAP20602] 
_i— Vipera berus CRISP [CAP74089] 



94r Opt 
97P Opt 

L-ri 



Ophiophagus hannah ophanin [AA062996] 
Ophiophagus hannah crisp-b (Indonesia) 
Naja kaouthia kaouthin-1 [ACH73167] 
Naja atra natrin-1 [Q7T1K6] 
53] Dispholidus typus CRISP-DIS2 [AAZ75596] 

Telescopus dhara CRISP-TEL1 [AAZ75606] 
' — Trimorphodon biscutatus CRISP-TRI1 [AAZ75607] 
Phiiodryas olfersii CRISP-PHI2 [AAZ75604] 
Rhabdophis tigrinus tigrin [Q8JGT9] 



Leioheterodon madagascariensis CRISP-LEI1 [AAZ75600] 
jEnhydris polylepis CRISP-ENH2 [AAZ75599] 
-s 1 Enhydris polylepis CRISP-ENH1 [AAZ75598] 

Eublepharis macularius crisp (SAL, SK) 

iAno//s carolinensis piscivorin-like [XP 003215375] 



100 'Anolis carolinensis crisp {testes) 
— Eublepharis macularius crisp (SAL) 
•Anolis carolinensis helothermine-like [XP 003215243] 
99r Heloderma suspectum CRISP [ACE95060] 
1 1 
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Heloderma horridum helothermine [AAC59730] 
40/ Varanus varius CRISP-VAR10 [AAZ75609] 
Varanus varius CRISP-VAR11 [AAZ75610] 
Varanus komodoensis CRISP [ABY89664] 

50i Varanus varius CRISP-VAR9 [AAZ75608] 

fif- Varanus acanthurus CRISP-VAR1 [AAZ75587] 
■ Varanus acanthurus CRISP-VAR2 [AAZ75588] 



40J | 

1 1 



90| — 



— Varanus acanthurus CRISP-VAR6 [AAZ75592] 
Varanus acanthurus CRISP-VAR4 [AAZ75590] 
Varanus acanthurus CRISP-VAR5 [AAZ75591] 
Varanus acanthurus CRISP-VAR8 [AAZ75594] 
Varanus acanthurus CRISP-VAR7 [AAZ75593] 
Varanus acanthurus CRISP-VAR3 [AAZ75589] 
Homo sapiens CRISP1 [NP 001192149] 



82 I 

65 1 



Mus musculus Crispl [NP 033768] 

Mus musculus Crisp3 [NP 033769] 

- Rattus norvegicus Crispl [NP 074050] 
Homo sapiens CRISP3 [EAX04348] 
Homo sapiens CRISP2 [AAI07708] 

1 — Mus musculus Crisp2 [NP 001191000] 
94I — Rattus norvegicus Crisp2 [NP 001011710] 



Supplementary figure S16 

Maximum likelihood tree of crisp sequences. Genbank accession numbers are given in square brackets. Sequences from this 
study are in orange text and the location of transcripts is provided in bracketed blue text (VG, venom gland; SAL, salivary 
gland; SK, skin; SCG, scent gland; AG, accessory gland; VMNO, vomeronasal organ). Numbers above branches 
correspond to Bootstrap values for 500 replicates. The * denotes the gene duplication event that produced crisp-a (shaded 
purple) and crisp-b (green). 
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Notechis scutatus venom C-type lectin [ABP941 19] 
Pseudechis porphyriacus venom C-type lectin [ABP941 08] 
Cryptophis nigrescens venom C-type lectin [ABP94123] 

■ Suta nigriceps C-type lectin [ACY68719] 
Pantherophis guttatus (SCG) 
Opheodiys aestivus (SAL, SCG) 
Thrasops jacksoni Lectoxin [A7X3Z0] 

Echis coloratus ctl-a (VG, SCG) 
Pantherophis guttatus (SAL, SCG) 
Bungarus fasciatus C-type lectin [Q90WI8] 

■ Pseudechis australis venom C-type lectin [ABP94084] 

MoreKa sp/fofa lectin-Mor-1 JAGI971 52] 
j — Pantherophis guttatus (SAL) 

1 — Micrurus corallinus C-type lectin precursor JACS74990] 
Oxyuranus scutellatus venom C-type lectin [ABP94096] 
Hoplocephalus stephensii venom C-type lectin [ABP94114] 
Bungarus flaviceps c-type lectin [ADF50042] 
Opheodrys aestivus (SAL) 
Leioheterodon madagascariensis Lectoxin [A7X401] 
Pantherophis guttatus (SAL, SCG) 
Opheodrys aestivus (SAL, SCG) 
Opheodrys aestivus (SAL, SCG) 
Pantherophis guttatus (SAL, SCG) 
Hoplocephalus stephensii venom C-type lectin [ABP94112] 
Oxyuranus scutellatus venom C-type lectin [ABP940931 
Pseudechis porphyriacus venom C-type lectin [ABP94128] 
Pseudechis australis venom C-type lectin [ABP94101] 
Oxyuranus microlepidotus venom C-type lectin [ABP94125] 

94 ■ Macrovipera lebetina C-type lectin B8JB4XT07] 

Macrovipera lebetina C-type lectin B6 [B4XT05] 
Echis coloratus ctl-b (VG) 

Deinagkistrodon acutus Agkicetin-C subunit beta [Q9DEA1] 

Echis coloratus ctl-c (VG) 

Bitis arietans CTL-1 [Q6X5T5] 
Macrovipera lebetina C-type lectin B9 [B4XT081 
Macrovipera lebetina C-type lectin B1 [B4XT00] 

Echis coloratus ctl-e (VG) ^ 
— Echis multisquamatus EMS1 6 subunit A TQ7T2Q1] 
Daboia russeilii P31 alpha subunit jADK22831] 
■ Bitis gabonica C-type lectin [Q6T7B51 
Macrovipera lebetina C-type lectin A9 [B4XSZ9] 
Macrovipera lebetina C-type lectin A2 [B4XSZ21 
Macrovipera lebetina C-type lectin A1 0 [B4XSY5J 
m ■ — Crotalus adamanteus C-type lectin JAEJ319671 
Gloydius halys Coagulation factor IXIQ9DG391 
Tnmeresurus flavovindis Coagulation factor IX[P23806] 
" ■oioratus ctl-f (VG) 
_ coloratus ctl-g (VG) 
Deinagkistrodon acutus Agkisacutacin [Q9IAM1] 
Echis coloratus ctl-h (VG) 
Macrovipera lebetina Coagulation factor X-activating enzyme [Q696W1 ] 
Echis coloratus ctl-i (VG) 
Echis coloratus ctl-j (VG) 
Macrovipera lebetina C-type lectin A1 5 [B4XSZ0] 





38 | 
|5J1 




Macrovipera lebetina C-type lectin A1 2 [B4XSY71 
Daboia russeilii dabocetin alpha subunit [ADK22821] 
Deinagkistrodon acutus Agglucetin subunit alpha [Q8JIV9] 
Crotalus horridus CHH : B [P81508] 



Deinagkistrodon acutus Agkicetin-C subunit alpha [Q9DEA2] 
aestivus (SAL) 

■ Anoiis carbiinensis c-type lectin BiL-like [XP 0032239231 

■ Anoiis carolinensis c-type lectin BFL-2-like [XP 003223920] 
■ Homo sapiens C-type lectin [AAG0051 4J 



; _1| Homo sap. ... ^ ..... . 

—I ' Meieagris gaiiopavo c-type lectindomain family 2 member D-like [XP 003214104] 

1 Gallus gallus C-type lectin [NP 001026644] 

j Rattus norvegicus C-type lectin domain family 4 member K [XP 002726471] 

n_i Homo sapiens C-type lectin domain family 4 member K [NP 056532] 

31 1 Mus muscuius C type lectin [CAC82936] 



0.2 



Supplementary figure S17 

Maximum likelihood tree of c-type lectin (ctl) sequences. The putative viper venom-specific ctl clade is shaded green. 
Genbank accession numbers or Ensembl gene IDs are given in square brackets. Sequences from this study are in orange text 
and the location of transcripts is provided in bracketed blue text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, 
scent gland). Numbers above branches correspond to Bootstrap values for 500 replicates 
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g^iEchis pyramidum acidic phospholipase A 2 [AAN77203] 
Bitis arietans P LA2-1 8 partial [AAX86638] 
Echis carinatus acidic phospholipase A 2 [AAN77202] 
Echis coloratus phospholipase A 2 [AAK49823] 

£c/7/'s coloratus PLA 2 IIA-e (VG) 

Echis carinatus phospholipase A 2 [AAP41217] 




Cerastes cerastes phospholipases A 2 [AC092622] 
^^^^^ Bitis gabonica phospholipase A2-1 [AAR06850] 



Vipera aspis ammodytin 12(A) [CAE47206] 
Vipera beras ammodytin 12 [AAN59985] 
99p Echis coloratus phospholipase A 2 [AAK49822] 
Echis coloratus PLA 2 IIA-d (VG) 
Daboia russellii phospholipase A2-I [AAP1 3805] 
Macrovipera lebetina phospholipase A 2 [CAR40186] 

Vipera aspis ammodytin 11(C) [CAE47168] 
Lachesis muta phospholipase A 2 [ADB77855] 
Agkistrodon piscivorus phospholipase A 2 [AAB25034] 
— Sistrurus catenatus phospholipase A 2 [ABY77925] 

Gloydius halys phospholipase A 2 [AAB71 848] 
^Calloselasmarhodostoma phospholipase PLA 2 [AAF03254] 
Crotalus oreganus PLA 2 2 [AEY82464] 

Crotalus atrox acidic phospholipase A 2 [AAL36974] 
Protobothrops elegans phospholipase A 2 [BAE72889] 
Bothrops aspermyotoxin Iva [ABY55159] 
Vipera am mo dytes ammodytin L [AAB20876] 
1 00i Echis ocellatus acidic phospholipase A 2 [AAN 77204] 
Bitis arietans PLA2-32 [AAX86637] 
Echis ocellatus phospholipase A 2 [CAQ72890] 
Echis coloratus PLA 2 IIA-c (VG) 

Echis coloratus PLA 2 IIA-b (SCG) 
Echis coloratus PLAo IIA-a (SCG) 
Pelodiscus sinensis PLA 2 [ENSPSIG00000005632] 
1 00 ^Meieagrisgaliopavo PLA 2 group MA [ENSMGAG00000005585] 

— Gallus gallus PLA 2 group MA [ENSGALG000000 18977] 
991 Homo sapiens PLA 2 group MA [ENSG00000 188257] 

58^^^— Mus musculus type II phospholipase A 2 [AAB06315] 
Rattus norvegicus PLA 2 group MA [EDL80890] 
Petromyzon marinusPLA 2 group MA [ENSPMAG00000005846] 



0.2 



Supplementary figure S18 

Maximum likelihood tree of Group IIA Phospholipase A 2 {PLA 2 Group IIA) sequences. Genbank accession numbers or 
Ensembl gene IDs are given in square brackets. Sequences from this study are in orange text and the location of transcripts 
is provided in bracketed blue text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, scent gland). Numbers above 
branches correspond to Bootstrap values for 500 replicates. The Echis coloratus sequences for PLA 2 IIA-a and IIA-b are 
likely to be alleles or due to individual variation between samples used in this study, and probably represent the same gene. 
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ar r Pyihan regius Serine protease (SCG) 

~~ '■Python rnotunis Serine protease {Pooled tissue) 
^j-DeinagUistrodon aetrftrs acubin [CAB4G431] 
IrLachesis stenophrys serine protea se[ABB7G2 79] 
jCratatus adamanteus Snake venom serine proteinase Ta [AFJ4B2G0] 
: ^Srstrcrrtrs eeferrattrs edwardsi serine proteinase isoform 5 [ABG2G971] 
-Profobotrrrops flavavfridrs Kallikrein [013053] 

jrlWacrovrpera febefrira serine protease VLSP-1 precursor [ADNQ491 S] 

Echis aoloraius serine protease [ADI47572] 

r jEah}spyrantsdum feaheyf serine protease [ADI47547] 
J ^cfris cannabis socbureki serine protease [ADI47550] 
'—Macrovipera rebetfna Kallikrein [Q5JH55] 
rEerrrs ocerrattrs serine protease [ADE4514Q] 

Echis oceJJatus sen ne protea se [ADE4513B] 
r-Echis cotoratus serine protease [ADI4757G] 

Cerastes cerastes Kallikrein [07SYF1] 

■i&jrrrs pyrami dum Jeatteyi serine protease [ADI47553] 
^cfrfs ocettaOis serine protease [ADI47556] 

irCrotatus adamanteus Snake venom serine proteinase 3b [AFJ49255] 
\Zrotatus dunssus terrifious gyroKin-like B1 7 serine protease precursor [ABYSSBS'I] 
-Gloydius ifi^tfrferrsrs thrombin -I ike eniymc ussurase [AAL4-5222] 
'^rotobotfrrops mucrosquamatus serpentokallikrein-'l [AAG27254] 
Crotalus adamanteus serine proteinase 5 [AEJ31 999] 
Zrotatus adamanteus Snake venom serine proteinase 1 [AFJ492 52] 
roydfifsbaryssalmobin [AACG133H] 

ijDeinagtfistrodon acutus snake venom serine protease Da-3G [AFR11 355] 
^■DerrrajjffrstrodDJT acutns thrombin-like protein 1 [AAVV5GG0B] 
■Vfrido'/fpera srejnecjerT Kallikrein [Q71 OH 6] 
■Botftrops atrox batroKobin [AAA43553] 
iomrops jararaca Kallikrein [Q9PTU3] 
=i— VTpera rrrfrorskff putative serine protease nikobin [CBW30775] 

Causus rhombeatus Kallikrein-Cau3 [ABU6B559] 

Betas coloratus Serine protease d (VG) 

itoydius shedaoensis Thrombin-like enzyme gloshedobin [PQC5B4] 
rtoydius ussuriensis serine protease gussurobin [Q3UVX1] 
■rotates adamanteus serine proteinase 9 [AEJ32003] 
- LBotfrrops jararacussu serine protease [AAP4241 6] 
■Viridovipera stejnegeriit&nom serine protease 1 [AAN5234G] 
■izCrotatus adamanteus Snake venom serine proteinase 5[AFJ49253] 
''VCrotalus atrox Kallikrein [Q3QHK2] 
j ' — Bomrops jararacussu serine protease [ABC24637] 
-Bottirops asperthrombin-like eniyme [ABB762S0] 

i 1 Dothrops jararaca platelet aggregating serine peptidase [AAQS25S0] 
'\aches;smuta Kallikrein [Q27J47] 



Echis cotoratus serine protease [ADI475GS] 
Crotalus adamanteus seri ne protei na se 7 [AEJ 32001] 
Crotatus adamanteus Snake venom serine proteinase G [AFJ49259] 
Gloydius nalys pa I la bi n [CAA04G 1 2] 
TGIoydius halys pa I la se [AAC 343 93] 
Gtoydius frafysagkiripin [AQK44E57] 
dr-Daboia russerrii srsrrrerrsrs RW-V gamma-like protein precursor [ADP3355B] 
iWacnsvrpera lebetina Factor V activator [Q9PT41] 
Echis cannabis soe Fri/rePf r seri ne protea se [ADI47575] 
Echis colorahis Seri ne protea se f |VG \ 
Echis coioratus serine protease [ADI47574] 
Macrwfpera lebemra Kalli krei n [Q9 PT40] 
Qitis gabonica Kallikrein [OGTG37] 
-Causus rfronrbeattrs Kallikrein-Cau1 [ABUG3557] 
- fEchis coioratus serine protease [ADI475G1] 
- \Ecnis oceflatus serine protease [ADE45141] 
Echis oceJlatns serine protease [ADI47564] 
-Gioydius brevrcaudus serine protease [AAD01 623] 
■ ijCrotatus atrox Serine protease catroxase-1 [QE3QHK3] 
ICrotalus adamanteus Kallikrein [Q98UUK2] 
-Frotobothrops mucrosquamatus Kallikrein [Q91 51 0] 
rOphiophagus barman- Kallikrein [EF0S0S37] 
* jLapemfs barduvfcffff serine protease precursor [AAV933G7] 
gBungarus muKcinctus Kallikrein [EF03033] 
tlaja atra Kallikrein [EF030334] 
>hi!odryas olfersff kalli krei n -Phi4 [AAZ75G2 3] 
Trfrodryas olfersrf serine protease precursor [ABI74G94] 
^h'rtodryas orfersir kallikrein-Phi3 [AAZ75G27] 
^frfrodryas orfersir kallikrein -Phi2 [AAZ7562S] 
"hiiodryas offersij kallikrein-Phil [AAZ75G25-] 
■PFrrrodryajr orrensFf ka I li krei n-Ph i 5 [AAZ75629] 
YTbamnophis eisgans Serine protease (Pooled tissue) 
HTranrrropFrfs sfsgans Serine protease {Pooled tissue) 
^Tiamnophis etegans Serine protease {Pooled tissue) 
^r-Fogoira vitticeps Serine protease {Brain) 

~l — Anotis carorfrrerrsfsttirombin-like enzyme stejnefibrase-1 -like [XPO 0 322 772 5] 
^^H/aranus indicus kalli krei n-Vind1 [ADK39243] 
^-TH/aramrs kDnrodoerrsfs: Kallikrein [BSCJU5] 
' — Varanus konrodoeirsrs: Kallikrein [B6CJUS] 
'-Varanus erenTfuskallikrein-Verel [ADK39263] 

Cirbfepri aris macut anus Serine protease (SAL, SCO, SK) 

■z:f\bronia ^ranrfrrea kallikrein-Abr-1 [AFU63206] 
"iGerrnanahis infematis kallikrein-Ginfl [ADK39259] 

- I Hefodenna suspecbim efrretuirr kallikrein toxin 1 [ACE950S4] 

1 -i— SerrfforroftrjB fnfenrarFjs kallikrein-Ginf3 [ADK392G7] 
'-Gerrfronotus firfenTaris kallikrein-Ginf2 [ADK392G0] 

McnQderpFrr:5: donrestfca cationic trypsin-like [XP001 3G2G77] 

Mus muscu!us tryps in 10 precu rs or [N PO 0 1 0 340 3 5] 
flafliTS rrorve^fcus anionic trypsin-2 precursor [NP036SS1] 
ftonro sapi ens Protease serinel {trypsin 1) [AAI23227] 
—Taenropygia guttata trypsin l-P1-like[XP0021 9501 3] 
Gaiiusgalius trypsin II-P29 precursor [TJP99071 5] 
'■Meteagris gattopavo tryps i n ll-P29-like [XPQ032Q2672] 

■f-sff rnerf a obatumnae Th rom bi n [EN S LAC G0000001SSOS] 

—Xenopus tropicaVss Thrombin [EN SXETG00000001 932] 
■^r-Taeniopygia guttata Th rom bi n [EN STG UGOOOQ00105HO] 
Hjif ereairrfs gaUopavo Th rom bi n [EN 3 MG AG 00000010 924] 
: . YSattusgallus Thrombin [ENSGALG00000003332] 
: -Ferodrsctrssrireirsfs Thrombin [EH SP3IG0000000500G] 
i-^Analis carQrfrrensfsThrom bi n [EN SAC AG 00000003353] 
t^Mus ntuscutus Th rom bin [EN 3 MU 5 G 0 0 0 0 0 027249] 
^PflaftiTS rrDrveofcirs Thrombin [ENSRNOG0000001 632S] 
'-Homo sapferrs Thrombin [EN SG000001 3021 0] 

Monoderprrfs domestj ca Thrombin [EN SMODG0000001 9393] 



Supplementary figure S19 

Maximum likelihood tree of serine protease (sp) genes. Genbank accession numbers or Ensembl gene IDs are given in 
square brackets. Sequences from this study are in orange text and the location of transcripts is provided in bracketed blue 
text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, scent gland). Numbers above branches correspond to 
Bootstrap values for 500 replicates. 
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mam yEchis ocettatus group III snake venom metallo proteinase [ADW54354] 

Echis ocettatus group III snake venom metalloproteinase [ADW54355] 
is coloratus metalloproteinase [ADI47G45] 

Echis coloratus svmp-a (VG) 
Echis coioratus metalloproteinase [ADI47G1 0] 

Echis ocettatus group III snake venom metalloproteinase [ADW54348] 
Echis coioratus metalloproteinase [ADI47625] 

Echis coloratus metalloproteinase [ADI4761 6] 
Daboia russsUU hemorrhagic metalloproteinase [AAZ39880] 

Echis ocettatus group III snake venom metalloproteinase [ADW54342] 

Echis coloratus metalloproteinase [ADI47654] 

MjEchis coloratus svm p-d {VG) 

^MEchis coioratus metalloproteinase [ADI47G58] 

Echis coloratus metalloproteinase [ADI47620] 




'Echis coioratus metalloproteinase [ADI47635] 

^hI'j ridovipera stejneueri stejnihagin-E [ABA4Q759] 
iTrimeresurus ftavoviridis hemorrhagic metalloproteinase HR1b [BAB92D14] 
i Bothrops jararaca hemorrhagic metalloproteinase HF3 [AAG48931] 
mCrotalus adamanteus Snake venom metalloproteinase (type III} 7 [AFJ49245] 
mBothrops jararaca bothropasin precursor [AAC61986] 
mCrotalus durissus metalloproteinase P-IIIIABA42117] 
mCrotatus adamanteus Snake venom metallo proteinase (type II) 6 [AFJ49228] 
mCrotalus adamanteus Snake venom metalloproteinase (type II) 8 [AFJ49230] 
mGloydius hatys metalloproteinase [AALG0587] 
bgkistrodon contortrix metalloproteinase VMP-I precursor [ACV8393D] 
lAgkistrodon contortrix metalloproteinase [AAC59703] 
mCrotatus durissus metalloproteinase P-ll [ABA42116] 
Qgkistrodon contortrix metalloproteinase [AAC597D4] 
mCrotalus molossus Snake venom metalloproteinase [G8JJ51] 
^gkistrodon piseivorus metalloproteinase VMP II [ACV8393G] 
■ Bomrops asper type I metalloproteinase [ABB76282] 
'Bothrops asper type II metalloproteinase [ABB 76281] 
mmmmmmEchis carinatus metalloproteinase [ADI47724] 
Echis coloratus metalloproteinase [ADI47686] 
Echis ocettatus Group I snake venom metalloproteinase [CAJ01686] 
Echis coloratus metalloproteinase [ADI47693] 

Echis pyramidum metalloproteinase [ADI477GG] 
Echis coloratus metalloproteinase JADI477011 

miEchis carinatus metalloproteinase [ADI47725] 
Echis ocettatus group II snake venom metalloproteinase [ADW54326] 
Echis pyramidum metalloproteinase [ADI477G1] 
■icoiorafussvmp-hfSCG) 

Echis coloratus metalloproteinase [ADI4770] 
Vi ridovipera slejnegeri venom metalloproteinase [ABC73D79] 

m Echis cotoratus metalloproteinase [ADI47653] 
chis coloratus svmp-k fVG) 
^j^^^Echfs coloratus svmp-l {SK) 
^^^^"Echfs cotoratus svmp m {VG) 
^^^^"Echfs coloratus svmp n (VG, 5CG) 
■l/jpera ammodytes H3 metalloproteinase [AGL45259] 
mEchis coloratus svmp-o (SCGJ 

mEchis coloratus svmp-p (VG) 




B mmmmEchis carinatus metalloproteinase [ADI47727] 

^Ech' 



hEchis cotoratus svmp-q (VG) 
V mmF Crotatus adamanteus Snake venom metalloproteinase (type III) 2c [AFJ49235] 
mmmmmm^^Crotatus adamanteus Snake venom metalloproteinase (type III) 8 [AFJ49246] 
I mmm A gkistrodon piseivorus metalloproteinase VMP III precursor [ACV83935] 

/mp-r{SCG^ 
lEchfs cotoratus metalloproteinase |ADI476Q8j| 

'SCG) 

m Bitis arietans metalloproteinase MP-2 [AAX86G34] 
m Cerberus rynchopsvenom metalloproteinase [ADJ51D55] 

Phitodryas otfersii snake venom metalloprotease [ACS74986] 

Ala/a atra metalloproteinase atrase B [ADG02948] 
■Ala/a mossambica mocarhagin 1 [AAM51550] 
Bungarus fasciatus metalloproteinase [ABN72536] 
iDemansia vestigiata metalloproteinase [ABK63559] 
Notechis scutatus scutatease-1 [ABQ01138] 
Pseudechis australis australease-1 [ABQ01134] 
Ophiophagus hannah metalloproteinase [ABM87941] 
Thamnophis elegans (Pooled tissue) 
Naja atra metalloproteinase [ABM72547] 

Anolis carotitiensis /ENSACAGQODD0011428] 
Eublepharis macularius {SAL, SCG] 
on regius {SCG} 
Echis cotoratus svmp -u (SCG] 

Thamnophis elegans (Pooled tissue) 
Opheodrys aestivus (SCG) 
Pantherophis guttatus (SCG) 
Gattus gattus ADAM28 [ENSGALG0OOOOGGQ357] 
Me/eagrisga»opai/oADAM28[EHSMGAG00000001087] 
Taeniopygia guttata ADAM 28 [EN STG UG 00000 QQ4505] 
Monodelphisdomestica ADAM2B [ENSMODGQQ0QQQ1G767] 
'Homo sapiens ADAM 28 [ENSG00000042980] 

mmmwRattus norvegicus ADAM28 [EH SRNGGGG00001451 8] 
iMus musculusADAM28 [ENSMUSGEKHJ0D0147251 
Monodetphis domestica ADAM7 [ENSMODG00000016755] 
Homo sapiens ADAM7 [ENSGOO00OOG920G] 

Rattus norvegicus ADAM7 [ENSRNOG00000014054] 
Mus musculusADAM7[ENSMUSGOD00002205Gl 
iMonode^hrs domesrica Decysin 1 [ENSMODG00000016762] 
-Homo sapiens Decysin 1 [ENSG0Q0001 34028] 

'Mus muscutus Decysin 1 [ENSMUSG00000022057] 
mRattus norvegicus Decysin 1 [ENSRNQGOOQ0003Q724] 




Supplementary figure S20 

Maximum likelihood tree of snake venom metalloproteinase (svmp) genes. Genbank accession numbers or Ensembl gene 
IDs are given in square brackets. Sequences from this study are in orange text and the location of transcripts is provided in 
bracketed blue text (VG, venom gland; SAL, salivary gland; SK, skin; SCG, scent gland). Numbers above branches 
correspond to Bootstrap values for 500 replicates. 

20 



10 20 30 40 50 60 70 

| | | | | | | | | | | | | | 

GD441513.il RNase-Gillf2 ATGGCTGCAA AGGGATCCTG TGTGTTCTGG TTCCTTTTGG CATCCGCCTG GATAGTGATG AAGAGTGATG 
GD441515.il RNase-Gillf4 



80 30 100 110 120 130 140 

| | | | | | | | | | | | | | 

GD441513.il RNase-Gillf2 CAGCTGACAC CTTTGAATCC T T CAAAGAGC TGCACGTTGA TTATCCAAAG ACTGGAGCTC CAAATGATGA 
GD441515.il RNase-Gillf4 



ISO 160 170 1B0 190 200 210 

| | | | | | | | | | | | | | 

GD441513.il RNase-Ginf2 TGAATACTGC AAAAAAGVVA VGGGGGGAAG AGGCCGGACC AAACTCAAGG CCAATACTTA CATTCATGCC 
GD441515.il RNase-Ginf4 



220 230 240 250 260 270 280 

| | | | | | | | | | | | | | 

GD441513.il RNase-Gillf2 CCAGACAGTG AACTTTTAGC CGCCTGCAAT CGGAAGAAGT ATAAACTTAA TCATGAATAC GGCAGGACCT 
GD441515.il RNase-Gillf4 



290 300 31D 320 330 340 350 

| | | | | | | | | | | | | | 

GD441513.il RNase-Gillf2 CTCGCCTTCC AACCACTTTG TGTACGTATG GTGATAGAGT ATTCTTGGOA AGCTCCTTGC CAGGAACG&A 
GD441515.il RHase-Gillf4 -- 



360 370 380 390 400 

| | | | | | | | | | | 

GD441513.il EHase-Gillf2 TAAGGTTCTG TGTGTGAATT GGAAGCCAGT TGCATTTAGA GGATTCAACG CCTAA 
GD441515.il EHase-Gillf4 



Supplementary figure S21. Alignments of Gerrhonotus infernalis ribonuclease sequences [3] 
GU441513 ("Gin£2") and GU441515 ("Ginf4") showing 100% similarity at the nucleotide level. 
Sequence lengths reflect those on Genbank and have not been edited. 
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10 20 30 40 EG 60 70 



GU441517.1 RNase-Gltlf6 ATGGCTGCAA AGGGATCCTG TGTGTTCTGG TTCCTTTTGG CATCCGCCTG GATAGTGATG AAGAGTGATG 
GU441516.1 RNase-Ginf5 



BO 90 100 110 120 130 140 



GU441517.1 RITase-Gitif 6 CAGCTGACAC CTTTGAATCC TTCAAAGAGC TGCACGTTGA TTATCCAAAG ACTGAAGCTC CAAATGATAA 
GD441516.1 RNase-Ginf5 



ISO 160 170 ISO 190 200 210 



GD441517.1 RKase-Cinf 6 TGAATACTGC AAAAAAGTTA TGGGGGGAAG AGGCCAGACC AAACTCAAGG CCAATACTTA CATTCATGCC 
GD441516.1 RNase-Gmf5 



220 230 240 230 260 270 2S0 



GD441517.1 ENase-Ginf6 CCAGACAGTG AACTTTTAGC CGCCTGCAAT CGGAAGAAGT ATAAACTTAA TCATGAATAC GGCAGGACCT 
G0441516.1 RNase-Ginf5 



290 300 310 320 330 340 350 



GD441517.1 Rtlase-Ginf6 CTCGCCTTCC AACCACTTTG TATGGCTGCA AAGGGAT CCT GTGTGTTCTG GTTCCTTTTG GCATCCGCCT 
GU441516.1 ENase-Gitif5 



3 SO 



GD441517.1 RNase-Gitlf6 GCATAGTGAT GAAGAGTGA 
GD441516.1 ENase-Ginf5 



Supplementary figure S22. Alignments of Gerrhonotus infernalis ribonuclease sequences [3] 
GU441516 ("Ginf5") and GU441517 ("Ginf6") showing 100% similarity at the nucleotide level. 
Sequence lengths reflect those on Genbank and have not been edited. 
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DQ139885 , 
DQ139884, 
DQ139886, 



DQ139885 . 
DQ139884, 
DQ139886 , 



DQ139885, 
DQ139884, 
DQ139B86 , 



DQ139885, 
DQ139884 , 
DQ139B86, 



DQ139B85, 
DQ139B84 . 
DQ139B86, 



DQ139885 , 
DQ139B84, 
DQ139B86, 



DQ139B85, 
DQ139884, 
DQ139B86 , 



DQ139B85, 
DQ139B84, 
DQ139886, 



DQ139B85, 
DQ139884 , 
DQ139886, 



1 CRISP-VAR5 
1 CRISP-VAR4 
1 CRISP-VAR6 



1 CRISP-VAR5 
1 CRISP-VAR4 
1 CRISP-VAR6 



1 CRISP-VAR5 
1 CRISP-VAR4 
1 CRISP-VAR6 



1 CRISP-VAR5 
1 CRISP-VAR4 
1 CRISP-VAR6 



1 CRISP-VAR5 
1 CRISP-VAR4 
1 CRISP-VAR6 



1 CRISP-VAR5 
1 CRISP-VAR4 
1 CRISP-VAR6 



1 CRISP-VAR5 
1 CRISP-VAR4 
1 CRISP-VAR6 



1 CRISP-VAR5 
1 CRISP-VAR4 
1 CRISP-VAR6 



1 CRISP-VAR5 
1 CRISP-VAR4 
1 CRISP-VAR6 



10 20 30 40 BO 60 70 

| | | | | I I I I I I I I I 

AT GAT CC T GC TCAAACTGTA TTTGACCCTA GCTGCAATCT TATGTCAATC CCOTGGCACG ACTTCTCTTG 

SO 90 100 110 120 130 140 

| | I I I I I I I I I I I I 

ATCTTGATGA TTTGATGACT ACCAACCCTG AGATACAAAA T GAG AT CAT T AACAAGCACA ATGACCTACG 

130 ISO 170 180 130 200 210 

| | | | | | | | | | | | | | 

GAGAACAGTG GATCCCCCAG CTAAAAACAT GCTGAAGATG TCCTGGGACA ACATCATTGC AGAGAGTGCC 

220 2S0 240 2i0 26D 270 2S0 

| | | | | | | | | | | | | | 

AAACGTGCAG CACTGAGATG CAACCAAAAT GAGCACACAC CTGTCTCGGG AAGAACAATA GGTGGTGTGG 

290 S00 310 320 330 340 350 

| | | | | | | | | | | | | | 

TGTGCGGAGA AAATTACTTC ATGTCAAGTA ACCCTCGCAC ATGGTCTTTC GGCATTCAGA GTTGGTTTGA 

3€0 370 3S0 390 400 410 420 

| | | | | I I I I I I I I I 

TGAAABGAAC TACTTTAAGT TTGGTTTTGG ACCAACAAGA GCAGGTGTCA TGGTTGGCCA TTATACCCAG 

430 440 450 460 470 4S0 490 

| | | | | | | | | | | | | | 

GTGGTCTGGT ATAAGTCTTA CAAAATGGGA TGTGCGATCA ACTTGTGCCC TAATGAGCCC CTGAAGTACT 

AA. .AAC.T. G. .G.C.GG. AGT . 

500 510 520 530 540 550 560 

| | | | | | | | | | | | | | 

TCCTGGTTTG CCAGTACTGC CCAGGAGGGA ACG77G7AGG CCGGAAGTAT GAACCCTATG CAATCGGAGA 

670 530 B90 600 610 620 630 

| | | | | | | | | | | | | | 

ACCATGTGCA GCTTGCCCCA ACAA77GTGA CAACGGACTG VGCACTAACC CC7GTGAGCA CAGCAATCAA 



DQ139885 , 
DQ139884, 
DQ139B86, 



1 CRISP-VAR5 
1 CRISP-VAR4 
1 CRISP-VAR6 



640 660 660 

| | | | | I 

TACATCAACT GCCCAGATTT AACAAAACAG 



Supplementary figure S23. Nucleotide alignments of Varanus acanthurus CRISP sequences [4] 
DQ139885 ("CRISP-VAR5"), DQ139884 ("CRISP-VAR4") and DQ139886 ("CRISP- VAR6") 
showing almost total similarity at the nucleotide level. Sequence lengths reflect those on Genbank and 
have not been edited. 
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DQ139B88 , 
DQ139887, 
DQ139883, 



CRISP-VAR8 
CRISP-VAR7 
CRISP-VAR3 



10 20 30 40 BO 60 70 

| | | | | I I I I I I I I I 

AT GAT CC T GC TCAAACTGTA TTTGACCC7A GC7GCAA7C7 7A7G7CAA7C CCG7GGCACG ACTTCTCTTG 



DQ133888 , 
DQ139887, 
DQ139B83, 



CRISP-VAR8 
CRISP-VAR7 
CRISP-VAR3 



SO 90 100 110 120 130 140 

| | | | | I I I I I I I I I 

ATCTTGATGA TTTGATGACT ACCAACCCTG AGATACAAAA TGAGATCATT AACAAGCACA ATGACCTACG 



DQ139888, 
DQ139887, 
DQ139883 , 



CRISP-VAR8 
CRISP-VAR7 
CRISP-VAR3 



150 ISO 170 180 130 200 210 

| | | | | | | | | | | | | | 

GAGAACGGTG GA7CCCCCAG CTAAAAACAT GC7GAAGA7G 7CC7GGGACA ACACCATTGC AGAGAGTGCC 



DQ139888, 
DQ139887, 
DQ139883 , 



CRISP-VAR8 
CRISP-VAR7 
CRISP-VAR3 



220 230 240 260 260 270 280 

| | | | | | | | | | | | | | 

AAACGTGCAG CACTGAGATG CAACCAAAA7 GAGCACACAC CTGTCTCGGG AAGAACAATA GGTGGTGTGG 



DQ139888, 
DQ139887 , 
DQ139883, 



CRISP-VAR8 
CRISP-VAR7 
CRISP-VAR3 



290 300 310 320 330 340 360 

| | | | | | | | | | | | | | 

TG7GCGGAG AAAATTACT TCATGTCAAG TAACCTTCGC ACATGGTCTT TCGGCATTCA GAGTTGGTTT 



DQ139888, 
DQ139887 , 
DQ139883, 



CRISP-VAR8 
CRISP-VAR7 
CRISP-VAR3 



360 370 3S0 390 400 410 420 

| | | | | | | | | | | | | | 

GATG&AAGGA ACTACTTTAA GTTTGGTTTT GGACCAACAA GAGCAGGTGT CATGGTTGGG CATTATACCC 



DQ139B88, 
DQ139887, 
DQ139883, 



CRISP-VAR8 
CRISP-VAR7 
CRISP-VAR3 



430 440 450 460 470 4S0 490 

| | | | | I I I I I I I I I 

AGGTGGT'CTG G7ATAAGTCT TACAAAATGG GATGTGCGAT CAACTTGTGC CCTAATGAGC CCCTGAAGTA 



DQ139888 , 
DQ139887, 
DQ139883, 



CRISP-VAR8 
CRISP-VAR7 
CRISP-VAR3 



500 510 520 530 540 550 560 

| | | | | I I I I I I I I I 

C77CC70O77 7GCCAG7AC7 GCCCAGGAGG GAACG77G7A GGGGGGAAG7 A7GAACGC7A 7GCAA7CGGA 



DQ139888 . 
DQ139887, 
DQ139883 , 



CRISP-VAR8 
CRISP-VAR7 
CRISP-VAR3 



570 580 530 600 

| | | | | | | | . . 

GAACCATGTG CAGC77GCCC CAACAACTGT GACAATGGAC TG 



Supplementary figure S24. Nucleotide alignments of Varanus acanthurus CRISP sequences [4] 
DQ139888 ("CRISP-VAR8"), DQ139887 ("CRISP-VAR7") and DQ139883 ("CRISP- VAR3") 
showing almost total similarity at the nucleotide level. Sequence lengths reflect those on Genbank and 
have not been edited. 
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10 20 30 40 50 €0 70 

| | | | | | | | | | | | | | 

GU441521.1 ESP-Vk.oml A7GGCATC77 "CA7AGC7"" CC77G7GTAC AGCTGGG77C 7CC7A7CGCT TGCAGCAGGT CAAGTTACCT 
EU1954 62.1 MMP2-Var2 

30 90 100 110 120 130 140 

| | | | | | | | | | | | | | 

GU441521.1 ESP-Vkoml ATCCAGAA7C 7CCAGCC7G7 AAG777CCCT TCATTTATGA GGGCAAAGCC TTCACTACAT GCACGGAATA 
EU1954 62.1 MMP2-Var2 

ISO 160 170 180 130 200 210 

| | | | | | | | | | | | | | 

GU441521.1 ESP-Vk.oml 7GG7"C7ACT GACAAAAC7C CA7GG'7G7GC CACAACCTCA AAC7A7GACA GGGATCGTAA A7GGAAGCCG 
EU1954 62.1 MMP2-Var2 

220 230 240 250 260 270 280 

| | | | | | | | | | | | | | 

GU441521.1 ESP-Vk.oml 7G7GC7G7CA AAGAG7A7GG AGAAC7GAGT GTCAACTACC CTGAAGATCC ATGTAGATTT CCTTTCACCT 
EU19S4 62.1 MMP2-Var2 

230 300 310 320 330 340 330 

| | | | | | | | | | | | | | 

GD441521.1 ESP-Vkoml ACAAAGGCAA GACTTACTCT GGTTGTACAG GAGCTQGGAG AGAGGATGGG AAGCTTTGGT GCTCTATCAG 
EU1954 62.1 MMP2-Var2 

360 370 3E0 330 400 410 420 

| | | | | | | | | | | | | | 

GD441521.1 ESP-Vk.oml CAAAAATCAT GATGACAACT CACAATTGGT ATTCTGTGAG CCTTCAGATC CAGCCCCCTG CTACTTTCCT 
EU1954 62.1 MMP2-Var2 

430 440 450 460 470 480 480 

| | | | | | | | | | | | | | 

GD441521.1 ESP-Vk.oml TTCAAATACA AGCAAAAATC CTACTCTGAC TGCACCATGA ATGGGAGTTT TGATGGGCAT TGGTGGTGTG 
ED1954 62.1 MMP2-Var2 

500 310 520 530 540 550 560 

| | | | | | | | | | | | | | 

GU441521.1 ESP-Vk.oml CCACAACCGC AGA77A7GAC AAGGACAGCA AGTGGAAGGC AT GT AT TACT GAAGAATATG GAGGAAACTC 
EU1954 62.1 MMP2-Var2 

670 580 530 600 610 620 630 

| | | | | | | | | | | | | | 

GD441521.1 ESP-Vtoml GAA7GGTGGG CAG7G7G7G7 7CCCC77CAC C7ACATGGAC AAGGAATATA CTACATGCAC CAA7GAGGAT 
ED1954 62.1 MMP2-Var2 

640 650 660 670 680 630 700 

| | | | | | | | | | | | | | 

GD441521.1 EEP-VtOMl ACAGA7GGGA GA77C7GG7G 7GCCACTACG GCAAACTATG ACAAGGACAA GAAA7GGAGC 77C7G7GC7G 
EU1954 62.1 tMP2-Var2 



GU441521.1 ESP-Vkoml 
ED1954 62.1 MMP2-Var2 



710 720 730 740 750 760 770 

| | | | | | | | | | | | | | 

ACACAAGAAT GAGCAAGGGG GAGATTGTCA AGCCAGGAAA TGTCGGCAGT TCTAAGGCGC CATGTTCCTT 



GU441521.1 ESP-Vk.oml 
ED1954 62.1 tMP2-Var2 



7G0 790 800 B10 820 830 840 

| | | | | | | | | | | | | | 

CCCCTTCATC TACAAAGGCA AGACTTACAC AGAGTGTACT T CCAAGGGGA CTTTAACATG GAAGCTCTGG 



GD441521.1 ESP-Vk.oml 
ED1954 62.1 >MP2-Var2 



850 860 870 8S0 890 900 910 

| | | | | | | | | | | | | | 

TGCTCTCTCA CCAGCAACTA TGATAGAGAT AAGAGATGGA GATACTGTCA GCCTTCAGAT TTAGAAACAC 



GU441521.1 ESP-VTioml 
EU1954 62.1 MMP2-Var2 



920 930 940 950 960 

| | | | | | | | | | | . 

AAAATCAGTT GCAGTCTAAT GAAGAGAAGC AAACAATAAA TATCAGCAAG GAATAA 



Supplementary figure S25. Alignments of Varanus komodoensis epididymal secretory protein 
(GU441521, "ESP-Vkoml") and matrix metalloproteinase (EU195462, MMP2-Var2) sequences [3] 
showing 100% similarity at the nucleotide level. Sequence lengths reflect those on Genbank and have 
not been edited. 
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10 20 SO 40 SO 60 70 



GD441522.1 ESP-Vk.om2 ATGGCATCTT TCATAGCTTT CCTTGTGTAC AGCTGGGTTC TCCTATCGCT TGCAGCAGGT CAAGTTACCT 
E0195454.1 MMP2 Varl 



30 90 100 110 120 130 140 



GD441522.1 ESP-Vliom2 ATCCAGAATC TCCAGCCTGT AAGTTTCCCT TCATTTATGA GGGCAAAGCC TTCACTACAT GCACGGAATA 
E01S5454.1 MMP2 Varl 



150 160 170 ISO 190 200 210 



GD441522.1 ESP-Vk.om2 TGGTTCTACT GACAAAACTC CATGGTGTGC CACAACCTCA AACTATGACG GGGATCGTAA ATGGAAGCCG 
E0195454.1 MMP2 Varl 



220 230 240 2S0 260 270 2S0 



GU441522.1 ESP-Vk.om2 TGTGCTGTCA AAGAGTATGG AGGTAATTCC AATGGAGCGC CATGCACCTT TCCCTTCATT TACCTTGGTC 
EU135454.1 MMP2 Varl 



290 300 310 320 330 340 350 



GD441522.1 ESP-Vk.om2 GTACTTACTA CACCTGCACC AACAAACTTG AGCATAAAGG ACGGTACTGG TGCGCCACAA CAGGGAGCTA 
ED195454.1 MMP2 Varl 



360 370 380 390 400 410 420 



GD441522.1 ESP-Vk.om2 TGATAAGGAC CGGAAGTGGA GTTTCTGTGC AG AC AT C AAA CTGAGTGTCA ACTACCCTGA AGATCCATGT 
ED135454.1 MMP2 Varl 



430 440 450 460 470 480 490 



G0441522.1 ESP-Vk.om2 AGATTTCCTT TCACCTACAA AGGCAAGACT TACTCTGGTT GTACAGGAGC TGGGAGAGAG GATGGGAAGC 
ED135454.1 MMP2 Varl 



500 510 52 0 530 540 550 560 



GU441522.1 ESP-Vk.om2 TTTGGTGCTC TAT CAGCAAA AATCATGATG ACAACTCACA ATTGGTATTC TGTGAGCCTT CAGATCCAGC 
ED135454.1 MMP2 Varl 



570 580 590 600 610 620 630 



GU441522.1 ESP-Vk.om2 CCCCTGCTAC TTTCCTTTCA AATACAAGCA AAAATCCTAC TCTGACTGCA CCATGAATGG GAGTTTTGAT 
ED195454.1 MMP2 Varl 



640 650 560 670 680 690 

| | | | | | | | | | | | 

G0441522.1 ESP-Vk.om2 GGGCATTTGT GGTGTGCCAC AACCGCAGAT TATGACAAGG ACAGCAAGTG GAAGGCATGT 
ED135454.1 MMP2 Varl T 



Supplementary figure S26. Alignments of Varanus komodoensis epididymal secretory protein 
(GU441522, "ESP-Vkom2") and matrix metalloproteinase (EU195454, MMP2 Varl) nucleotide 
sequences [3] showing near total sequence identity (lbp difference over 675bp of aligned sequence). 
Sequence lengths reflect those on Genbank and have not been edited. 
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